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Cosmological magnetic fields induce temperature and polarization fluctuations in the cosmic mi¬ 
crowave background (CMB) radiation. A cosmological magnetic field with current amplitude of 
order 10~ 9 G is detectable via observations of CMB anisotropies. This magnetic field (with or with¬ 
out helicity) generates vector perturbations through vortical motions of the primordial plasma. This 
paper shows that magnetic field helicity induces parity-odd cross correlations between CMB temper¬ 
ature and B-polarization fluctuations and between E- and B-polarization fluctuations, correlations 
which are zero for fields with no helicity (or for any parity-invariant source). Helical fields also con¬ 
tribute to parity-even temperature and polarization anisotropies, cancelling part of the contribution 
from the symmetric component of the magnetic field. We give analytic approximations for all CMB 
temperature and polarization anisotropy vector power spectra due to helical magnetic fields. These 
power spectra offer a method for detecting cosmological helical magnetic fields, particularly when 
combined with Faraday rotation measurements which are insensitive to helicity. 

PACS numbers: 98.70.Vc, 98.80.-k, 98-62.En 


I. INTRODUCTION 

and isotropic 
the Universe 
On the other 

hand, the presence a cosmological magnetic field |2|, Lil Li 1^1 IZt ti| alters this situation. In this paper we focus on 
vector perturbations induced by a cosmological magnetic field, existing since radiation-matter equality or earlier (for 
a specific inflation model see Q)- Both an homogeneous magnetic field, as well as the more realistic stochastic one, 
induce transverse MHD (Alfven) waves jj[. Non-decaying cosmological Alfven waves result in CMB temperature and 
polarization anisotropies 0, B, fi. I7L IlflUllj. 3 

A cosmological magnetic field could have helicity M- Magnetic field helicity plays an important role in the MHD 
dynamo, used in some models of galactic magnetic field amplification llffl . Unlike gravitational waves which are 
damped on scales below the Hubble radius at decoupling (the corresponding damping multipole number l ~ 100), 
Alfven waves survive down to smaller damping scales (l ~ 2000) fzj, thus a cosmological magnetic field can affect small- 
scale CMB fluctuations. From this point of view, the vector mode is more relevant for constraining a cosmological 
magnetic field from CMB fluctuation measurements. 

In this paper we present analytic expressions for all CMB fluctuation vector power spectra that arise from a 
helical cosmological magnetic field. In particular, we also compute how magnetic helicity affects parity-even CMB 
fluctuation power spectra. We propose a scheme to constrain cosmological magnetic helicity from CMB temperature 
and polarization anisotropy observations. The symmetric part of the magnetic held spectrum can be reconstructed 
from measurements of the rotation of the CMB polarization plane as a consequence of the Faraday effect m ; this is 
because magnetic helicity does not contribute to the Faraday rotation effect I18lll9ll2fl . On the other hand, the helical 
part of the magnetic held spectrum induces parity-odd cross correlations between temperature and B-polarization 
anisotropies, and between E- and R-polarization anisotropies ESQ; such cross correlations are not induced by the 


The linear theory for the evolution of small inhomogeneities in the standard spatially homogeneous 
cosmological model, developed by Lifshitz and others 'l|, shows that vector perturbations 1 decay as 
expands. 2 This is why vector perturbations are usually neglected when computing CMB fluctuations. 


* Electronic address: tinatin@phys.ksu.edu 

'Electronic address: ratra@phys.ksu.edu 

1 Vector perturbations are also known as transverse peculiar velocity or vorticity perturbations. 

2 This is for adiabatic initial conditions, which are consistent with observations. However, a non-decaying vector perturbation exists if 
there is an initial large-scale photon-baryon fluid vorticity compensated by a neutrino vorticity such that the total large-scale vorticity 
in relativistic components vanishes 0. Such a mode can be excited only after neutrino decoupling 0- 

3 In addition to the vect or mode, a cosmological magnetic field generates CMB fluctuations through the scalar and tensor modes as 

described in Refs. IE E3 [3 [l3 • 
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Faraday effect f20j . 4 

For our computations we use the formalism of Ref. [g|, extending it to account for magnetic field helicity. To 
compute CMB temperature and polarization anisotropy power spectra we use the total angular momentum method 
of Ref. |23|. Our results are obtained using analytic approximations and are valid for l < 500. We present results in 
terms of a ratio between CMB fluctuation contributions from the symmetric and helical parts of the magnetic field 
power spectrum. 

The rest of the paper is organized as follows. In the next section we model the helical magnetic field source term for 
vorticity (Alfven) waves. In Sections III and IV we present analytic approximations for the parity-even and parity-odd 
CMB fluctuation power spectra contributions induced by magnetic helicity. In Section V we discuss our results and 
conclude. In Appendixes A and B we present some details of the computation including analytic approximations used 
for some of the integrations. 


II. MAGNETIC FIELD INDUCED VECTOR PERTURBATIONS 


A. Magnetic field source term 


We assume the existence of a cosmological magnetic field generated during or prior to the radiation-dominated 
epoch, with the energy density of the field a first-order perturbation to the standard Friedmann-Lemaitre-Robertson- 
Walker homogeneous cosmological spacetime model. Neglecting fluid back-reaction onto the magnetic field, the 
spatial and temporal dependence of the field separates, B(t, x) = B(x)/a 2 ; here a is the cosmological scale factor. As 
a phenomenological normalization of the magnetic field, we smooth the field on a comoving length A with a Gaussian 
smoothing kernel oc exp[— a; 2 /A 2 ] to obtain the smoothed magnetic field with avera ge v alue of squared magnetic field 
B\ 2 = (B(x)-B(x))|a and magnetic helicity If a 2 = A|(B(x)-[VxB(x)])|a. See Ref. |2(| for a more detailed discussion. 

We also assume that the primordial plasma is a perfect conductor on all scales larger than the Silk damping 
wavelength A s (the thickness of the last scattering surface) set by photon and neutrino diffusion. We model magnetic 
field damping by an ultraviolet cut-off wavenumber fco = 2tt/Xb 1T3 . 


kn 


Mpc 
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k\ 
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riB+ 3 


h. 


(i) 


Here ns is the spectral index of the symmetric part of the magnetic field power spectrum (see Eq. 0 below), h is 
the Hubble constant in units of 100 km sec -1 Mpc -1 , k\ = 2n/X is the smoothing wavenumber, and Xd A s- This 
assumes that magnetic field damping is due to the damping of Alfven waves from photon viscosity. 

Assuming that the stochastic magnetic field is Gaussianly distributed, and accounting for the possible helicity of 
the field, the magnetic field spectrum in wavenumber space is m, 


(B*(k)Bj(k')} = (2n) 3 6^(k - k')[P„(k )P B {k) + teyifciPff (*:)]• (2) 


Here Py(k) = Sij — kikj is the transverse plane projector with unit wavenumber components ki = ki/k, etji is the 
antisymmetric tensor, and (k — k') is the Dirac delta function. We use 


BA k) 


d 3 xe lk ' x Ijj(x), 


P?( x ) 


d 3 k 
(27r) 3 


BA k), 


( 3 ) 


when Fourier transforming between real and wavenumber spaces; we assume flat spatial hypersurfaces. Ps(fc) and 
Pn{k) are the symmetric and helical parts of the magnetic field power spectrum, assumed to be simple power laws 
on large scales, 


9—2\3o 2 o 2\3rr2 

and vanishing on small scales when k > kp. Here T is the Euler Gamma function. These power spectra are generically 
constrained by Ps(fc) > \Pn(k)\ mm/ which implies nn > tib 00]. In addition, finiteness of the magnetic field 


4 This is not true for a homogeneous magnetic field l21l . 
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energy density requires n B > —3 (to prevent an infrared divergence of magnetic field energy density). Finiteness of the 
magnetic field average helicity requires n B > —4; this is automatically satisfied as a consequence of Tin > n B > —3. 

To obtain the magnetic field source term in the transverse peculiar velocity perturbation equation of motion we 
need to extract the transverse vector part of the magnetic field stress-energy tensor Ty(k). This is done through 
Iljj(k) = (Pjh(k)fcj + Pjb (k) hi ) k a Tab (k) , and the Ily tensor is related to the vector (divergenceless and transverse) 
part of the Lorentz force L^(k) = kjHij (k) = P ib (k)k a T a b(k) (Eq. (2.16) of Ref. j|). For the normalized Lorentz 
force vector IR = l\ V ^ /k, the general spectrum in wavenumber space is similar to Eq. ©; that is, 5 


(n*(k)IT, (k')} = (27r) 3 <5 (3) (k - k') 


Pijf(k)+ie ijq k q g(k) , 


( 5 ) 


where f(k ) and g{k ) represent the symmetric and helical parts of the vector Lorentz force power spectrum. These 
functions are related to the magnetic field stress-energy tensor spectrum through, 


(27r) 3 <^(k - = -[P aj (k)^P am (k')fc( + P Ql (k)%Pa/(k')C]( t* (k)r im (k')), 

(27r) 3 <5 (3) (k - k')g(k) = --k q [ejmqkik'i + ejz g fcjP(J(T*-(k)T/ m (k')). 


The functions f(k) and g{k) are evaluated in Appendix A, 


f(k) ~ k D ) 2nB + 3 
g(k) ~ g\k{Xk D ) nB+nH+2 


1 + 


n B 


ub + 3 \koJ 
1 + ""- 1 


k \ 2nB+3 


n B + 3 \kn 


- T h {X k D ) 2nH+3 

ns+riH+2' 


1 + 


n H 


- 1 / k \ 2nn+3 


n H + 4 \k D J 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


for k < fcu; f(k) and g{k) vanish for k > k B - Here Ph and Q are constants that depend on the magnetic field 
power spectrum indexes n B and , and the power spectrum normalization, 


XFb 


__ 

16(2n B + 3)r 2 (n B /2 + 3/2)’ 


XFh 


24(2n B + 3)r 2 (n ff /2 + 2)’ 


__ 

24(n B + n H + 2)T(n B /2 + 3/2)T(n H /2 + 2)' 


( 10 ) 


The contribution of magnetic field helicity to the symmetric source f(k) is negative (as in the case for tensor pertur¬ 
bations (3)- The magnetic source terms in Eqs. m and © vanish on scales smaller than the cutoff scale A < X B 
because of magnetic field damping. The singularities at n B = —3/2, n B = —3/2 and n B + n B = —2 in Eqs. (ITHll are 
removable F° r n B > —3/2 the terms proportional to k^ B+3 and k^ H+3 dominate in the expression for f(k) 

in Eq. ©, and the symmetric source term depends on the cutoff wavenumber k B but not on k , and so is a white 
noise spectrum [23j . 


B. Vorticity perturbations 

A cosmological magnetic field contributes, via the linearized Einstein equations, to all three kinds of perturbations, 
scalar, vector, and tensor modes (for a recent review see 0 )- Here we focus on the effects a stochastic magnetic 
field with helicity has on vector perturbations. 6 The vector metric perturbation may be described in terms of two 
gauge-invariant divergenceless three-dimensional vector fields, the vector potential V (which is a vector perturbation 
of the extrinsic curvature), and a vector parameterizing the transverse peculiar velocity of the plasma, the vorticity 
O = v — V, where v is the spatial part of the four-velocity perturbation of a stationary fluid element [3j. In the 
absence of a source a vector perturbation decays with time and so can be ignored. 


5 We thank A. Lewis (private communication 2004 and Ref. Q) for pointing out a missing (27r) 3 factor in the expression for (II*(k)IIj(k')) 
given in Ref. ;6||. In what follows we use other results from Ref. §| corrected for a similar missing factor. 

6 See Ref. O for a study of helical vorticity fields. They did not account for magnetic field helicity acting as a source in the helical 
vorticity field perturbation equation of motion. 
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Since electromagnetism is conformally invariant it is possible to rescale fields by appropriate powers of the scale 
factor and simply obtain Maxwell’s equations in the expanding universe from the Minkowski spacetime Maxwell 
equations. Since the fluid velocity is small the displacement current in Ampere’s law may be neglected; this implies 
the current J = V xB/(47t). The residual ionization is large enough to ensure that magnetic field lines are frozen into 
the plasma, so the induction law takes the form (9/i9t)B = V x (v x B). As a result the baryon Euler equation for 
v has a Lorentz force L(x) ~ — (B(x) x [V x B(x)]} /( 47 t) as a source term. The photons are neutral so the photon 
Euler equation does not have a Lorentz force source term. The Euler equations for photons and baryons are p, ^, 


U 7 + f(v, 


7 


flh 


a r> T ( 

- 5 (v ’ 


Vb) 

Vb) 


= 0, 

Il v >(k) 

a 4 (pb +Pb)’ 


( 11 ) 

( 12 ) 


where an overdot represents a derivative with respect to conformal time 77 , and = v 7 — V and fib = v;, — V are 
the vorticities of the photon and baryon fluids. Here f = n e aTd is the differential optical depth, with n e the free 
electron density and ax the Thomson cross section, R = (pb +Pb)/(p-y + p 7 ) — 3pb/4p 7 is the momentum density 
ratio between baryons and photons, and ^ is the transverse vector (divergenceless) part of the Lorentz force. 

The average Lorentz force (L(x)) = — (B x [V x B])/(47t) vanishes, while the r.m.s. Lorentz force (L(x) • L(x)) 1 / 2 is 
non-zero and acts as a source in the vector perturbation equation. The magnetic helicity spectrum Pff(fc) contributes 
to the symmetric part of the Lorentz force spectrum, see the expression for f(k ) in Eq. |5|l. As a result magnetic 
helicity affects the symmetric vorticity perturbation spectrum via the Euler equation, and so contributes to parity- 
even CMB fluctuations. The helical part of the Lorentz force spectrum is completely determined by g{k), Eq. fU|l. 
and acts as a source for the helical part of the vorticity perturbation spectrum. Solving the Euler equations in the 
tight-coupling limit when v 7 ~ v&, we have 0I1, 7 


H (± 1 ) (^,k) 


m( ±1 )(k)?7 

(1 + i?)(p 7 0 +P70) 


(13) 


Here p 7 0 and p 7 o are the photon pressure and energy density today, and the component vanishes due to the 
transversality condition. This result can also be obtained from the Einstein equation jf|. 

The expression in Eq. G3J is valid on scales A larger than the comoving Silk scale A s- On smaller scales the photon 
viscosity becomes comparable to the effects of the magnetic field and so must be accounted for in the Euler equation. 
On these smaller scales with k > ks 0, 


tt (± 1 ) (??,k)~ 


n (±1 >(k) 


(fcL 7 /5)(p 7 o + p 7 o) ’ 

where L 7 is the photon mean free path length. 

We define the CMB anisotropies in terms of the vorticity perturbation power spectrum , 8 

<L!*(k)%(k')} = (27 t) 3 <5 ( - 3 )(k — k') [P y |n| 2 (fc)+*fc,ey,w(fc) 


(14) 


(15) 


Here |fi| 2 (fc) is the symmetric part of the vorticity power spectrum, directly related to the symmetric magnetic field 
source spectrum f(k) = |n| 2 (fc), while u>(i],k) is the helical part of the vorticity spectrum and is determined by the 
helical magnetic field source term g(k), 


ki 7 


f2(77,fc)~J (1 + R)(pjo + Pjo) 


(fcL 7 / 5 )(p 7 o +P70) 


n(fc), k < k s , 

n(fc), k> ks, 


k ) ~ 


krj 


1 2 


(1 + R) (p 7 o + P 70 ) 

1 

(fcL 7 / 5 )(p 7 o +P70). 


g{k), k < k s , 

2 (16) 

g(k), k > k s . 


Here the factor 1 + R for k < ks reflects the suppression of the vorticity field due to the tight coupling between 
baryons and photons, because photons, being neutral, are not influenced by the magnetic field Lorentz force source. 


7 We use the “helicity” basis of Sec. 1.1.3 of Ref. 153 and decompose a vector A = 1 e ( 7 ), where 1 ^ and e® are the unit 

basis vectors. The unit vector is chosen to be in the direction of wave propagation e-T = k and e- ^ 1 ^ (k) = — i(e 1 ± ie 2)/\/2. 

8 Both the vorticity ft and the vector potential V appear in the equations for the CMB temperature and polarization fluctuations. The set 
of the equations governing the vector perturbation dynamics is given in Eqs. (35)-(38) of Ref. jjj- These show that the contribution from 
n to the CMB anisotropies dominates over that from V, justifing neglect of the V contribution to CMB temperature and polarization 
anisotropies HIH 
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III. PARITY-EVEN CMB FLUCTUATIONS FROM MAGNETIC HELICITY 


Cosmological magnetic field vector and tensor mode contributions to large angular scale (/ < 100) CMB fluctuations 
are of the same order of mag nitude while small angular scale (l > 100) CMB fluctuations are dominated by the 
vector mode contribution [14( . Stochastic non-helical magnetic field effects on the CMB temperature and polarization 
anisotropies are discussed in detail in Refs. Si0- Here we compute the CMB vector mode fluctuation arising from 
the helical part of the magnetic field power spectrum (also see Ref. M). 

To compute the CMB temperature and polarization anisotropy power spectra we use the total angular momen¬ 
tum representation |22[. Given CMB temperature and polarization anisotropy integral solutions of the Boltzmann 
temperature equation, the CMB fluctuation power spectra are [22}, 


cr '=- 


„ ,2 V - ' k) 

^ 2Z +1 21 + 1 


(17) 


where r/o is the conformal time now and X is either 0, E , or B , which represent respectively the temperature, E- 
polarization, and R-polarization anisotropies. For the vector mode the sum includes only terms with m = ±1. In 
what follows we use results from Ref. @ for the CMB fluctuation power spectra C x s x induced by the symmetric 
(non-helical) part of magnetic field power spectrum and proportional to /s(fc) ~ f d 3 p Pb(p)Pb (|k — p|). 

The complete parity-even CMB fluctuation power spectra may be expressed as, 


C? x ' 


— r '■ 

— G 


xx' 
( S)i 


-C 


XX' 
(A)l - 


(18) 


where the are the power spectra induced by magnetic helicity, i.e., proportional to fu(k). The minus sign reflects 

the negative contribution of magnetic helicity to the total parity-even CMB fluctuation power spectra. 9 This result 
also holds for the tensor mode case [f4} . The fractional difference n xx = 1 — C x ^ /C x s x , where 0 < n xx < 1, can 
be used to characterize the reduction of the ^parity-even CMB fluctuation power spectra amplitudes as a consequence of 
non-zero magnetic helicity. The ratio Cfy* /C x s x may be expressed in terms of g(k)/f(k), i.e., in terms of Pqh/Pqb 
and spectral indexes uh and ns- We find that for any parity-even CMB fluctuation power spectrum 


2(2 n B + 3) f PHok n D H ~ nB 
3(2 uh + 3) V Pbo 


2 

P-a X { n B, riH , l) ■ 


(19) 


Here 1Z XX (ubtI^HjI) are dimensionless functions of l and spectral indexes ns and nu, and the index a = 1 and 
a = 2 corresponds to uh > —3/2 and nu < —3/2, respectively. In this section we present explicit forms for 1Z XX . 
Parity-odd CMB fluctuation power spectra, such as C/ es and Cj ^ B , receive a contribution only from the helical 
part of the magnetic field source power spectrum, g(k ), i.e., from terms proportional to J d?pP b(p)Ph[\^ ~ p|) and 
/ d 3 pP h(p)Pb{\^- - p|). 


A. CMB temperature anisotropies 


Vector perturbations induce CMB temperature anisotropies via the Doppler and integrated Sachs-Wolfe effects. 
Neglecting a possible dipole contribution from the velocity perturbation v today, and using the Boltzmann temperature 
transport equation solutions 0; for vector perturbations (see Eqs. (5.2), (5.5), and (5.6) of Ref. 0), we get 10 


riO& 
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(kDV o) 


2njj+3 
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n H + 4 V k D 


/ \ 2n H +3l 

f J Hi/2(^o>- (20) 


9 Reference E3 ignores the magnetic helicity contribution to the symmetric vorticity power spectrum. 

ignore all terms proportional to f d^p Pjj (p)Ph (|k — p|). 

10 For C®,® see Eq. (5.7) of Ref. |]. 


That is, in their computation they 
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Here J;+ 1 / 2 ( 2 ;) is a Bessel function, and we have used the analogy with Alfven velocity, va\ = B\/ 47r(p 7 o +p 7 o), 
to introduce the “helicity” velocity vj/a = H\/ ^/4:Tt(p i0 + p 7 o)- For ti# > —3/2 the integral expressing the CMB 
temperature anisotropy in Eq. ca is dominated by the first term (1) in the curly brackets, while for nn < —3/2 the 
second term oc ( uh — 1 )/{nn + 4) in the curly brackets dominates. 

The integral in Eq. (®J) is split into two parts. The first integral is evaluated using the solution for vorticity that 
does not account for the effects of viscosity, Eq. d, while the second integral makes use of the solution in Eq. o 
which accounts for the effects of viscosity. However, numerical integration shows that the first integral J^ s dominates, 
since the second integral with vorticity damped by photon viscosity (when ks < k < fcp) contributes less than of order 
1 % for l < 500, for all values of the spectral indexes ub and nn 0- Thus the contribution of the second integral 
may be safely neglected. The first integral is evaluated using Eqs. (5.9)-(5.11) and (5.15) of Ref. JfJ, and Appendix 
B of this paper. 

Retaining only the leading-order terms for l < 500, for uh > —3/2 we find, 


nf 6 


1 , 

2 (jib + 3)(nB + 2 ) / \ 2 " B+3 

< n B \ks) 

2 2nB+e (riB + 3)F(—2n B — 4) ( k S fio\ / k D r, 0 \ 2nB+3 
n B T 2 (-n B - 2 ) \ l ) \ l ) 


ns > —3/2, 
—2 < n B < -3/2, 

—3 < ns < —2. 


( 21 ) 


Even for a magnetic field with maximal helicity (when |Ptr(fc)| = Pb(/c)), for ub — riH > —3/2 the ratio /C^®® = 
2/3 (/«p e = 1/3) and the CMB temperature anisotropy power spectrum C® e > 0. For ns > —2 the function 
'Bf e (nB,riH 1 l) is independent of Z. 11 For —3 < ns < —2, IZf e (ns, nn, 0 is a growing function of l, scaling as 
^-2n B -4 The maximum growth rate of l 2 occurs for ng = —3. 

For uh < —3/2 we distinguish two different regions, — 2 < nn < — 3/2 and — 3 < uh < — 2. Using Eqs. iBB m 
(also see Eqs. (5.10)-(5.12) and (5.15) of Ref. 6J), and b we find, 


7ep e 


(n g -l)(n a +3)(nB + 2) fk D \ 2{nB uh) 
n B (n H + i)(n H + 2) \ks J 

2 2nB+6 {n H - 1 ){n B + 3)r(-2n B - 4) /fc£\ 2(nfl " nH) f k s y o\ 
n B (n H + 4)(n ff + 2)r 2 (-n s - 2) \k s J \ 1 / 

(■ n H - 1 )(n B + 3)r(-2n_B - 4)F 2 (-n g - 2) / k D r] 0 


‘2^‘ n b—^uh — 1 

77, d (m tt 


2n B +4 


-2 <n B < -3/2, 
n B < -2 <n H < -3/2, (22) 


2(n B -n H ) 


The function IZf 0 is independent of l for —2 < ub < tih < —3/2, while grows as l -2 " 8-4 for —3 < ns < 
—2 < uh < —3/2, which coincides with the growth rate of lZf e (n,B,nH,l) for —3 < ub < —2 and uh > —3/2. 
Thus if ub < —2, independent of whether uh > —3/2 or —2 < nn < —3/2, lZ® e (n,B, nn, 0 oc ;- 2ra s- 4 . For 
—3 < ns < uh < —2 the function TZf & (nB, nu, l) monotonically increases with l if nn > ns, while it is independent 
of l for nn = ns- 

The approximate expressions above are accurate to better than 15 % for uh > —2, to better than 30 % for 
—2.5 < nn —2, and to better than a few percent for nn < —2.5 0. We do not reproduce the explicit forms of the 
power spectra for the cases nu = —3/2 or uh = —2, and ub = —3/2 or ns = —2, since the corresponding results can 
be easily derived via a straightforward extension of the computation presented in Ref. jjj; see Eq. (5.16) there. 


B. CMB polarization anisotropies 

Vector perturbations generate both E and B CMB polarization anisotropies |2^, [26|. Since scalar perturbations do 
not induce a magnetic ( B ) CMB polarization anisotropy, the future detection of a B polarization signal will indicate 
the presence of a vector and/or a tensor perturbation mode. In this subsection we consider the CMB polarization 
anisotropies that result from vector perturbations induced by an helical cosmological magnetic field. 


11 For ng = —2 there is a weak dependence on /. 7?.®®(ns, ng , l) ~ 1 /1 n(A', 9 ryri//)■ 
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1. E -polarization 


Using the integral solution for the CMB electric (E) polarization Ei , see Eq. (6.6) of Ref. |jj, the CMB .E-polarization 
power spectrum from magnetic helicity is, 12 

C EE = n 1UM ox (2tt ) 2nn+8 vjj X (kpr] 0 ) 2n »+ 3 2 ( 7?dec?70 V 

{A)l 1 )[ + , M(2n H + 3)T 2 (n H /2 + 2) (k xVo )^+ 5 ^ dec {l + R de J 

(23) 


fks 

/ dkk 5 

l ! n H - 1 f k N 

2jih+3 

) 

[(i + 

Ji+3/2(kilo)' 

'o 

bh + 4 \k D y 


[ K (kr]o) 2 

kilo 


The contribution from the integral (i.e., the contribution from the region where vorticity is damped by photons) 
is negligible. When evaluating the integral in Eq. 03 we retain only the dominant term, and this differs depending on 
whether n h > —3/2 or nn < —3/2. We also use Eq. IIB3I) to approximate the cross term Ji+i/ 2 (kr) 0 )Ji +3 / 2 (kr) 0 ). To 
evaluate TZ EE (ns, riH, l ) we consider the following three regions: (i) ub > —3/2 and nn > —3/2; (ii) —3 < tib < —3/2 
and n,H > —3/2; and, (iii) —3 < hb < —3/2 and nn < —3/2. 

For tih > —3/2, using Eqs. IB2I1 IB5I) and Eq. IB8I) and retaining only the leading terms in the limit where 

ksVo > h we find, 


Hf E 


1 , 

2 {tib + 3) 2 
n B 



ns > —3/2 
-3 <n B < -3/2. 


(24) 


The function 1Z EE = lZf e for hb > —3/2, and Hf E = (n,B+3)lZf e for —2 <hb < —3/2. The expression for lZf E is 
the same for all values of Ub in the range —3 < ns < —3/2, while TZf 0 differs depending on whether —2 < ub < —3/2 
or —3 < ub < —2, see Eq. (12111 . 

For tih < —3/2 we need to consider the ranges —2 < nn < —3/2 and nn < —2 separately 0. Using Eqs. (IB2I) IB5D 
and Eq. (□51, we find that for uh < —2 the expression for C E ^ L is identical to the one in the region —2 < nn < —3/2 
to within 15 % accuracy. This simplifies the computation, and for —3 < ub < nn < —3/2 we find, 


K EE (n B ,n H ,l) 


(n g -l)(n B + 3) 2 (ksY {nH - nB) 
n B {n H + 4)(n ff + 3) \k D ) 


(25) 


while for —2 < < uh < —3/2 we have lZ EE (nB, nn, l) = (ns + 3 ){n H + 2)/{{ns + 2 )(nn + 3))lZf e (ns, nn, l), 

where lZf e is given in Eq. (PI). 


2. B -polarization 


To compute the contribution from cosmological magnetic helicity to the CMB E-polarization power spectrum we 
use the integral solution for Bi , Eq. (6.16) of Ref. The CMB E-polarization power spectrum is, 

r BB (1 | n) (2? t) 2uh+8 V 4 hx (k D 7 lo ) 2nH + 3 2 ( ?7dec??0 X* 

{A)l [ ,[+) 54(2n H + 3)T^n H /2 + 2)(k x r, 0 )^H+e"f^{i + Rde J 


fks 

/ dkk 5 

\l+ nH ~] 

( k \ 2n ^+ 3 ' 

lo 

nn + 4 

{k D J 


J ;+i/ 2 ( fc? 7o) 

(ki 7o) 2 


(26) 


Here again the contribution from the region where viscous effects are important is negligibly small. An expression for 
C E S y is given in Eqs. (6.18)-(6.19) of Ref. 0. 

Using Eq. lB2ll and retaining the leading terms for nn < —3/2, we get, 


nf B ~ nf E , 


K% b ~ K% e , 


(27) 


12 Just as for the temperature integral solution (-)/. Ej is also expressed in terms of f2. To compute C^ E we use Eqs. im and CSJ along 
with Eq. 181 . An expression for is given in Eqs. (6.7) and (6.9)-(6.12) of Ref. |j|. 




















to better than 20 % accuracy 

Both functions R EE and R BB (i.e., the coefficients nf E and nf B ) are independent of l. This means that cosmological 
magnetic helicity reduces the CMB polarization power spectrum amplitudes by the same scale factor for the E- and 
the E-polarizations. The ratio between contributions to the E- or E-polarization signal from the helical and the 
symmetric parts of the magnetic field are independent of l for the entire range of spectral indexes in the case of the 
vector mode, while for the tensor mode case the ratios depend on l for uh < —2 [hi . 


3. Temperature-E -polarization cross correlation 

We may obtain the CMB temperature-E-polarization cross-correlation power spectrum Cf E from the integral 
solutions for the temperature and E-polarization anisotropies. Like C®® this power spectrum is also parity-even 
and only the symmetric part of the magnetic field source /s(fc) (which also contains a contribution from the helical 
part of the magnetic field power spectrum Pn(k)) contributes to it. As discussed in Ref. ‘22j (see Fig. 5 there), 
the vector dipole temperature anisotropy radial function j^ 1V ^ = ji^/l(l + \)/2/x does not correlate well with its 
.E-polarization anisotropy radial function = [ji/x 2 + j(/x\^/(l — 1 )(l + 2)/2 (here j[(x) is the partial derivative 
with respect to x of the Bessel function of fractional order jj), while the vector quadrupole temperature anisotropy 
radial function jj 2V ^ = \/3 1(1 + l)/2(ji/ x)' does. To compute the vector mode CMB temperature-E-polarization 
cross-correlation power spectrum we therefore have to retain the term proportional to jj 2V> in the vector temperature 
integral solution 0; which was neglected previously 13 in the derivation of C®® (see Eq. lt20l) l. With the j[ 2V ^ term 
the CMB temperature anisotropy integral solution is , 


Qi V \r,o,k) ~ 11(1 + 1 ) 


2 / + 1 


"^(?7dec5 &) 


jl(krj Q ) kL 7 

dec 


krjo 


(*-l) 


ji(km) ji+i(km) V 


(km) 2 


km j 


(28) 


is given in Eq. (7.3) of Ref. Jj, and for C^j l we get, 


Cffy = - VKl-W + W + Z) 


(27r) 


2n H +S „,4 
U H\ 


(kDV o) 


2uh+3 


18(2n ff + 3)T 2 (n ff /2 + 2) ( kxm ) 2nH+6 Vl + ^d, 


^7dec^?0 


Jf y dec 


r k s 


dk fc 4 



2rif/+3 


(i + 1 ) 




^+i/2 (km) Ji+i/2(km)Ji+3/2(km) 

(km) 3 (km) 2 

Ji + i/ 2 (kvo)Ji +3 / 2 (kvo) J 2 +3/2 (km) 


(km) 3 


(km) 2 


(29) 


The first two terms in the second pair of square brackets in this integral arise from the correlation of jj 1 ^ with . 
A numerical evaluation of the integral @ shows that these two terms roughly cancel each other as a consequence of 


the low correlation between j\ and e\ v |23|. This may also be seen by using Eqs. IB(il) and m, 

—— J l+i/2( x )- J l+i/2(x)Ji+3/2(x) ^ 7 }[Ji-i/ 2 (x)Ji+i/ 2 (x) - Ji+i/2(x)Ji+3/2(x)} ^ sin(2:r - Ztt) - sin(2x - hr - tt) ~ 0. 

(30) 

In what follows we neglect these two terms. We use Eqs. (IB2I) (IB5I) to evaluate the three terms in the curly brackets 
of Eq. . The teims from the correlation between j j and c j are suppressed by an additional factor of k dec , 
relative to the two terms from the correlation between jp V> and ep \ In the limit l 3> 1, these terms and the squared 
sum of the two Bessel function terms in the expression for C EE in Eq. 1231) (the last factor inside the integral of 
this equation) are almost identical. Thus apart from an overall minus sign, the C ® Fj are approximately equal to the 
corresponding C EE |g]. Here our approximation might not be as accurate because, accounting for the suppression 
factor &L 7 deo the neglected contribution from the correlation between and ep could be comparable to the 


retained contribution from the correlation between jl‘ 2V> and e 


(V) 


13 Since it is suppressed relative to the term proportional to j) 1 . 
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IV. PARITY-ODD CMB FLUCTUATIONS FROM MAGNETIC HELICITY 


Magnetic helicity induces parity-odd cross correlations between the E- and R-polarization anisotropies, as well as 
between temperature and R-polarization anisotropies mm- Such off-diagonal parity-odd cross correlations also 
occur in the case of an homogeneous magnetic field from the Faraday rotation effect fuj], but not in the case of 
a stochastic magnetic field, even one with non-zero helicity |2Cj. Faraday rotation measurements cannot be used to 
detect magnetic helicity Uil L1IJ, |2C|. A possible way of detecting magnetic helicity directly from CMB fluctuation data 
is to detect the above parity-odd CMB correlations or to detect the effects magnetic helicity has on parity-even CMB 
fluctuations. In this Section we study the parity-odd CMB cross-correlations generated from vorticity perturbations. 
The corresponding tensor mode contributions are derived in Sec. VI of Ref. |14ll . 


A. Temperature—B-polarization cross correlation 


To compute the cross correlation between the CMB temperature and E-polarization anisotropies we use the integral 
solutions for 0; and Bi given in Eqs. (5.6) and (6.16) of Ref. ji] and find, 


Cf B 


-y/(l-l)l(l + l)(l + 2) 


(27r 


\nB+riH+8 ,,2 


V A\ V HX 


(&£>%) 


n.B+HH+2 


18(n s + n H + 2)T(n B /2 + 3/2)T(n H /2 + 2) (, k x Vo ) nB+nH+6 


^dec^O 

1 + Rdet 


1 2 


L 


7 dec 


rks 


dk k 3 


Vo 


1 + 


n H ~ 1 
ns + 3 


k 


nB+riH+ 2 ' 


Ji 


1 + 1/2 


(ki ?o) • 


(31) 


To evaluate this integral it is helpful to consider separately the index ranges (i) n B + n B > —2 and (ii) uh+tib < —2. 

When n B + n B > —2 the integral on the r.h.s. of Eq. ED is dominated by the first term in the square brackets. 
Using Eq. (lB2l with p = 3 we find, 

^(ksvof ■ (32) 

Vo 

When riH + n B < —2 the integral on the r.h.s. of Eq. ED is dominated by the second term (oc k nB+na+2 ) in the 
square brackets. For ub+b-h < —5 the integral in Eq. ed </ 0 fe dkk nB+nH+5 J 2 +1 , 2 (kr]o)) converges for large ks which 
can then be extended to oo, and so the integral may be evaluated using Eq. (IBII) . We find for —6 < n B + tih < —5, 


c es = 


(27r 


\nB+riH+i bi* 


V A\ V H A 


{k D v o) 




27 (n B + riH + 2)r(ns/2 + 3/2)r(n^/2 + 2) ( k\po ) nB+nH+6 (1 + Rdec)vo 


Vde 


c? B 


in,B+riH+7 


(2 7r ) ,l B+ rl H'+ 8 2 7lB +"B+4 


A\ U H A 


9 (ns + uh + 2)T(n B /2 + 3/2 )T(tih/2 + 2) ( k\ijo ) nB+riH+6 


Vde 


1 2 


(1 + Rdec)V0 


L 


7 dec 


Vo 


n H — 1 \ r {-n B - n H - 5) 


ns +3 / F 2 (— ub/2 — uh/2 — 2) 


(33) 


For n B + > —5 (but still n B + < —2) the integral in Eq. (TIH) diverges at large ks and so the upper limit 

cannot be replaced by oo, and the integral cannot be evaluated by using Eq. ED- Instead we approximate it by 
using Eq. We find for — 5 < ns + uh < —2, 


Cf B 


(270"*+"*+ 7 VAX v 2 hx 


( ksvo) 


UB+nn+5 


9 (n B +n H + 2){n B + n H + 5)T(n B /2 + 3/2)T(n H /2 + 2) ( k\Vo) nB+nH+ 6 U 1 + Rdec)vo_ 


Vde 


- / 7 dec 

Vo 


\ tib +3 ) 


(34) 


When ns + tih = —5 the integration can be done by using Eq. lB2ll with p = 0. 

At large angular scales (l < 100) where the contribution from the tensor mode is significant, for n B + n h > —2 
the vector mode Cf B ^ X * and the tensor mode Cf B ^ (see Eq. (98) of Ref. [l^ l have the same l dependence oc l 2 . 
For all other values of spectral indexes n B and njj, the growth rate (with /) of Cf B ^ is faster than Cf B< - T \ In 
particular, when the integral in Eq. ED converges at large ks (for —6 < tib+Uh < —5), cf B ' yV ' > /Cf B ^ oc l 3 . When 
—5 < ns + tih < — 2, the integral for the tensor mode temperature-R-polarization cross-correlation power spectrum 
converges at large k (see Eq. (97) of Ref. j^|), while it diverges for vorticity perturbations, Eq. ED. resulting in 
C@b(v) / qQB(t) ^ j-nB-nff- 2 _ ra ^j 0 between temperature-B-polarization signals from vector and tensor modes 

is independent of the amplitudes of the average magnetic field (B\) and average magnetic helicity {H\). 
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For maximally helical magnetic fields with nn — ns, due to the suppression factor L 7i <j ec/vo the temperature- E- 
polarization cross-correlation power spectrum is smaller the temperature-B-polarization cross-correlation power 
spectrum C ® B , 14 Cf E <C C® B , but both are oc l 2 , if ns + uh > —5. The same suppression factor makes Cf 5 
smaller than C) ee . For an arbitrary helical field Cf* B /Cf E depends on the ratio (Pho/PB o)k'^ 1 ~ nB and order unity 
prefactors that depend on n B and A dependence on l appears only if ns + nn < —5 when the ratio C® B /Cf E 
decreases as oc l n B+n H +5^ 


B. E- and .B-polarization cross correlation 


To compute cross correlations between E- and B-polarization anisotropies we use the integral solutions for Ei and 
Bi given in Eqs. (6.14) and (6.16) of Ref. |fj. We find, 


cr 


-(1-W + 2 ) 


(27r) n B+ n H+ 


( k D rio ) 


UB+UH+ 2 


54(n s + n H + 2)r(n B /2 + 3/2)r(n B /2 + 2) ( k x r]o ) nB+nH+6 


'Hdec'HO 


1 + Rde 


T 2 

- L/ 7dec 


rks 

/ dk k 5 


( k \ nB+nH+2 ' 


/o 

ub + 3 

V k n ) 



(i + 1) 


^1+ 1/2 (^^o) Jl+l/2(k' , lo)Jl+3/2(krio) 

(A770) 3 (fc?7o) 2 


(35) 


The combination of Bessel functions in the second set of square brackets in this integral is identical to the first two 
terms in the second set of square brackets in Eq. 12611 , which is negligibly small according to Eq. 1301) . Also the E- and 
B-polarization anisotropy cross-correlation power spectrum has an additional suppression factor of fcL 7>c j ec relative to 
the expression in Eq. (1231) . This implies C EB <C Cf B . Note that this is consistent with the result of Ref. [§2j that 

j\ n does not correlate well with e[ V \ The corresponding C EB amplitudes in the tensor mode case (for l < 100) are 
not suppressed and are of the same order of magnitude as the tensor mode temperature-B-polarization anisotropy 
cross-correlation power spectrum 0- 


V. CONCLUSION 

In this paper we consider how cosmological magnetic helicity affects CMB fluctuations. Even for a cosmological 
magnetic field with maximal helicity, such effects may be detectable only if the current magnetic field amplitude is at 
least 10 -10 or ICC 9 G on Mpc scales. Our analytical expressions for CMB fluctuation power spectra are valid (to the 
accuracy of our approximations) for n B > —3. 15 

A cosmological magnetic field generates a .B-polarization signal via induced vector and/or tensor modes, so a 
detection of such a signal may indicate the presence of a cosmological magnetic field. 16 However, it has to be 
emphasized that a B-polarization anisotropy signal can also arise in other ways, such as from primordial tensor 
perturbations [25| , gravitational lensing |2f|, or Faraday rotation of the CMB anisotropy polarization plane 013 0 . 
The B-polarization anisotropy power spectrum l 2 C BB peak position may help identify the B-polarization source. For 
example, cosmological-magnetic-field-induced tensor perturbations only contribute on large angular scales l < 100, 
while B-polarization anisotropy from gravitational lensing has a peak amplitude l 2 Cp B ~ 10 -14 at l ~ 1000 [29[. The 
Faraday rotation B-polarization anisotropy signal from a field with B\ = 10~ 9 G (at A = 1 Mpc) and spectral index 
ub = —2 peaks at a substantially smaller scale l ~ 10 4 with a frequency-dependent peak amplitude l 2 C BB ~ 10~ 12 
(at 10 GHz) and l 2 C BB ~ 10 -14 (at 30 GHz) (20f. A non-helical cosmological magnetic field with B\ = 10 -9 G 
at A = 1 Mpc induces a B-polarization anisotropy signal via the vector perturbation mode with a peak amplitude 
l 2 C BB ~ lO" 13 at l 1000 |7]. We have shown that a magnetic field with maximal helicity results in the reduction of 


14 If n H ~n B > -3/2, C? E /Cf> B ~ L 7 ,dec/^o) 

15 How ever, a cosmological magnetic field with spectral index ub — 2 (as might be generated by an MHD cascade in the early universe 
m ), has significant power on small (galaxy cluster) scales and s o m easurements of Faraday rotation in clusters imply an upper limit 
on the smoothed amplitude B\ < 10 —12 G on Mpc scales mu Such a “blue” cosmological magnetic field cannot significantly 
affect CMB anisotropies. A strong constraint on magnetic field amplitude on Mpc scales for ub > —2 arises from gravitational waves 
production via a magnetic source, if the magnetic field is generated with a power law spectrum Ps(/c) = (k/k m ax.) nB Pb (kmax.) (where 
femax = ?7gen and r)g en is the moment of magnetic field generation) that can be extrapolated, without damping, down to the Hubble 
radius when the magnetic field is generated, e.g., « 10 — 4 M pc i f the magnetic field is generated at the elecroweak phase transition (or 
even smaller for a magnetic field generated during inflation) l28| . 

16 See Refs. l26l for CMB polarization anisotropy measurements. 
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the R-polarization anisotropy signal on all scales by a factor of 1/3 for —3/2 < ns — n h , relative to the non-helical 
magnetic field case. 

We have presented analytical expressions for all CMB fluctuation power spectra affected by cosmological magnetic 
helicity. Our results show that cosmological magnetic helicity can affect CMB anisotropies, in addition to the effects 
it has on MHD dynamo amplification and processes in the early universe HEHij- It would be useful to incorporate 
our analytical expressions for l 2 C xx into a numerical code (e.g., that of Lewis ;?]) to compute CMB temperature 
and polarization anisotropies generated by a general cosmological magnetic source. 

To set an observational limit on cosmological magnetic helicity H\ (or Pn{k)) one may proceed as follows. The 
first step is to determine the average magnetic field B\ (or Pg(fc)) using measurements of the Faraday rotation of the 
CMB polarization plane Il7l 11911201 . Then one may use measurements of the parity-odd temperature-R-polarization 
anisotropies cross-correlation power spectrum C® B for l > 100 (to insure that the tensor mode does not contribute). 17 
According to Eq. 423 . if B\ and ns are known, C® B is determined by H\ and nn (he., the helical part of the magnetic 
field spectrum). A future detection of 0 — B cross correlations may be used to constrain H\. It should be emphasized 
that on scales l > 100 magnetic-field-induced cross correlations between E- and S-polarization are negligibly small, 
C BB <C Cf B , which can be used as a cross-check of the source of B-polarization anisotropy. To bound the range of nn 
for a given tib, the l dependence of the ratio Cf^/Cq®® can be used. In particular, if, for —3 < ns < —2, C/®n; /^fs)i 

is an increasing function of l growing as l~ 2n B - 4 ; then tih > —2, while if it scales as l 2 ( n H~ n B) then ns < nn < —2. 
If, for tib < —2, C(^ 4 ®/C(®® is l independent, then uh — ns- If Ub > —2, C^®/C(®® is l independent for any allowed 
nn A tib > — 2. 

It is possible that there are other ways to detect magnetic helicity. On the other hand, Ref. 0 argues that a 
detection of magnetic helicity, even for cluster magnetic fields, is a very difficult task. 
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APPENDIX A: DERIVATION OF MAGNETIC FIELD SOURCE TERMS FOR VORTICITY 

PERTURBATIONS 


The derivation of the symmetric magnetic field source term for the vorticity perturbation equation of motion is 
given in jjj] . To obtain the complete magnetic field source term for the vector metric perturbation we use the Fourier 
transformed magnetic field energy-momentum tensor, 


,(k) = 


2(2t 


<Bp 


Bi(p)Bj(k - p) - ]-B l (p)B l (k - p)(5, 


(Al) 


and its wavenumber space power spectrum (r*- (k)-n m (k')) jf|. The symmetric part of the magnetic source power 
spectrum f(k) and the helical part g{k) of the magnetic source power spectrum can be obtained via Eqs. © and 0 , 
see Eqs. (A4)-(A7) of |j and Eq. (Al) of Ref. M- The part of f(k) = fB{k) — fnik) that depends on the symmetric 
part of the magnetic field power spectrum, /s(fc) ~ J d 3 p Pb(p)Pb( |k — p|), is jg, 


fs{k) 


16(2n s + 3)r 2 (n s +3/2) 


(A k D ) 2nB+3 + UB (A k) 2nB+3 
tib + 3 


(A2) 


The magnetic helicity contribution, fn(k) ~ f d 3 pP h( p)P h (\k — p|), to f(k) is, 


Mk) = 


1 

8(2tt) 5 


d 3 p P h (p)Ph (|k - p|) 


P( 1 -7 2 ) 


\/ k 2 — 2fcp7 + p 2 


(A3) 


17 Magnetic helicity is the only known source of diagonal temperature-5-polarization anisotropies cross correlations. A homogeneous 
magnetic field induces off-diagonal 0 — 5 cross cor relations E3, but not diagonal correlations. Faraday rotation does not induce 
diagonal or off-diagonal 0 — 5 cross-correlations EJ. 
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where 7 = k • p. The function f B {k) is positive but it contributes to f(k) with a minus sign, f(k) = f B Ck) — fn(k), 
and so it decreases the overall symmetric magnetic field source term. This result contradicts that of Ref. [jjj] . Reality 
requires P B (k) > \Pii{k)\, and the total symmetric source term is always positive, f(k) > 0. The mixed terms 
J d 3 p P B (p)P H (\k — p|) and f (Pp Ph(p)Pb( |k —p|) do not contribute to f(k), while the helical spectrum g{k) is 
completely defined by these terms, 


9{k) 


1 

16(27t) 5 


d 3 p P b (p)Ph{ |k 


(fc-2p7)(l-7 2 ) 

\Jk 2 — 2kpr/ + p 2 


(A4) 


The helical spectrum g(k) does not receive a contribution from the diagonal terms f d 3 p P B (p)P B ( |k —p|) and 

/ d 3 p P h {p)Ph{ |k - p|). 

To evaluate the expressions for f B {k) and g{k ), we first integrate over 7 and then integrate over p. The integration 
over 7 can be done analytically, see the Appendix of Ref. | 6 |- To integrate over p we approximate the result of 
the 7 integration by using the binomial expansion (1 + x) n = 1 + nx + n(n — l)x 2 /2 + 0(x 3 ). For the vector case 
the dominant terms for x <C 1 are those of quadratic order. Additionally, the integration is split into two parts, 
/ Q fc ° dp = / Q fc dp + J^ D dp, and the binomial expansion for k > p used for dp, while the second integral f£ D dp is 
evaluated using the binomial expansion for k < p ft| ij3| . The result is the vector perturbation source power spectrum 
symmetric ( f(k )) and helical ( g{k )) terms in Eqs. (JSJ) and © above. 


APPENDIX B: BESSEL FUNCTIONS INTEGRALS 


We need to evaluate integrals of the form f 0 ' s dx J p (ax)J q (ax)x~ b , which contain products of Bessel functions. For 
b > 0 when the integral converges and is dominated by x xs, the upper limit xs can be replaced by 00 (with an 
accuracy of a few percent for b > 1 , and 15-30 % for 0 < b < 1, depending on the value of p — q). We can then use 
Eq. (6.574.2) of Ref. 


/•OO 

/ dx J p (ax)J q (ax)x~ b = 

Jo 


a b - 1 T(b)T((p + q-b+l)/2) 


2 b T({-p + q + b+ l)/2)r((p + q + b + 1)/2)F((p - q + b + l)/2) ’ 


(Bl) 


which is valid for Re (p + q + 1) > Re b > 0, and a > 0. 

To evaluate the integral f^ s dx x p J 2 +1 ^ 2 (x) with p > 0 and xs l, we use the asymptotic expansion of J p [x) for 
large arguments, Eq. (9.2.1) of Ref. iu|, Ji+ 1 / 2 ( 2 ) ~ ^/2/(ttx) cos[a; — (l + 1 W/ 2 ] ~ y/2/{-Kx) cos [ 2 ; — (l + 1 )7 t/2 ]. 
Replacing the oscillatory function cos 2 by its r.m.s. value of 1/2, we obtain 'ft(l4l|. 


p > 0 

ln(xs/{l + 1 / 2 )), p = 0 . 


» J dxx"J^ 1/2 (x) ~ir^ 
irxJi +1 / 2 {x)Ji +3/2 (x) ~ 2 cos (x~(l + 1 )^j cos (x - (l + 2)^ =2 cos (x - (l + l)^j sin (x - (l + 1)^ 


(B2) 


For the integral J^ s dxx p Ji+\/ 2 {x)Ji +3 / 2 (x) we also use the large argument (x l) approximation for the Bessel 
functions, and find (see Ref. for a numerical check), 


= sin (2x — (l + 1 )tt) = (—l ) z+1 sin( 2 x). 


(B3) 


So for p > 0, 

r x s 


7r 


r x s r x s 1—lV / \ 

J dxx p Ji + i/ 2 (x)Ji +3 / 2 (x)~(—l) l+1 Ji dxx p ~ l sin ( 2 a;) ~—-—(a/g -1 sin( 2 a;z}) — (l + l/ 2) p_1 sin( 2 Z)J . 


(B4) 

This approximation tends to underestimate; it is good to a few percent for p > 1 and is within 30 % for 0 < p < 1. 
For p = 0 the integral j^ s dx Ji + \/ 2 (x) Ji +3 / 2 {x) may be evaluated using Eq. (11.4.42) of Ref. [31), 


rxs 

/ dx Ji+i/ 2 (x)Ji +3 / 2 (x) = 
Jo 


(B5) 


For the integral f u ' s dxx p [(l + l)Ji + i/ 2 {x)/x — Ji+ 3 / 2 {x )] 2 appearing in Eq. lt^3l) . for large enough l we have the 
approximation, 


(—^Ji+i/ 2 {x) - Ji+ 3 / 2 (x) S j « \J[+i/ 2 { x )} 2 — - Ji +3 / 2 {x)] 2 . 


(B 6 ) 
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We now approximate the cross term Ji-\/ 2 {x)Ji+ 3 / 2 {x) in the limit x^> l by using, 

TrxJi- 1 / 2 (x)Ji +3 / 2 (x) ~ 2cos (a; - l^j cos (z - (l + 2)^) = -2cos 2 (x - l^j . 

As in the computation of Eq. (IB2I we may replace the cos 2 by 1/2, and so find for p > 0, 

f xs fl + 1 \ 2 7r r x s / \ 

7T J dx X P y—^—Ji + i/ 2 {x) — Ji +3 / 2 (x) J — —J dxX P yJ l _ 1 i 2 — ‘2Jl-\/2Jl+3/2 J rJl+Z/2) 


. p t2 _ J X P o/p, p > 0 

dxX Jl “I ln(* s /(i + l/2)), p = 0. 


(B7) 


(B8) 
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